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■ Abstract Stationary differential systems with polynomial right sides are considered. 



Necessary and sufficient conditions are formulated when a given domain is a domain of 
asymptotic stability and the origin of coordinates is either focus or center. The problem 
of construction of a stabilizing control in the form of polynomial is studied. 
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^ 1 1 Introduction 

O ' 
(N , 

In the 19 th century the great French mathematician Henri Poincare formulated the problem 
^ \ of finding stability conditions for differential systems without calculation being a solution. 

At the end of 19 th and at the beginning of 20 th centuries the great Russian mathematician 
A.M.Lyapunov developed the mathematical stabilization theory for differential systems. 
For this purpose he developed two methods. The first one is based on the characteristic 
numbers of the fundamental matrix. The second one is based on construction of special 
functions called Lyapunov's functions which have properties similar to distance from a 
considered point to the origin of coordinates (it is supposed that the zero solution is 
stationary). Lyapunov's methods received world recognition. 

In this paper, differential systems with a right polynomial part /(.) of 

x = (x 1 ,x 2 , ...,x n ) G K n 

are considered i.e. 

f(x) = (f 1 (x),f 2 (x),...J n (x))\ 

where 

f (~\ — \ ^ „(?) Jl J2 _i n 

Jp\ x ) — a h,l 2 ,-,ln 1 2 ••• X n ! 
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(p) 

p G 1 : n, li, lii ■■■In are non-negative integers and * is the transposition sign, ct^ 2 ; n are 
real-valued numbers, I p is the set of degrees of the polynomial f p (x). 
To analyze the system 

x= f{x) 

a method is suggested which is different from Lyapunov's methods and based on a system 
transformation idea, so that we are able to say something definite about stability. 

From the technical point of view the important problem is finding a domain of asymp- 
totic stability. Conditions are given so that any domain including the origin of coordinates 
is one of them. Necessary and sufficient conditions are found when the origin of coordinates 
is either focus or center. 

(p) 

The case is considered when the coefficients a) ,, , depend on t. Sufficient conditions 
are given for the problems formulated above. 

The systems considered are interesting because asymptotic stability of the zero solution 
is equivalent to the following statement, that any solution of our system starting from any 
point of some small region of the origin of coordinates tended to be zero; this is not correct 
in general case. 

Further, we solve the problem of constructing a stabilizing control in any given region 
of the origin of coordinates for the system 

x= f(x,u), 

where x £ K n is a phase vector , u £ R r is a control, f(x, u) = (/i(x, u), fi(x, u), f n (x, u))* 
is a vector polynomial of x and u i.e. 



f p (x,u)= 

h,l2,.-,ln,mi,m2,—,m r elf 



a h,l2,...,ln,mi,m2,...,mr X l X 2 •" x n "l «2 — u r 



(p) 

p G 1 : n, l 1 ,l2,---,ln,m 1 ,m 2 ,...,m r are non-negative integers, a/i j / 2) ... ) j„ )jnijm2j ... jmr are 
real- valued numbers, I p is the set of degrees of the polynomial f p (x,u). We assume that 
the zero vector = (0, 0, 0) € 3i n is a solution of the system i.e. /(0, 0) = 0. 

The theorem was proved that any region including the origin of coordinates can be 
made a region of asymptotic stability if we choose a suitable control u(x). The control u{x) 
is chosen as a polynomial with a degree not higher than the degree of the vector-function 
f(x,u) as a function of x. 

The formulated problem about asymptotic stability of a system with a right polynomial 
part is very important for practical applications in physics and technics and can be found 
in the book by V.I.Zubov "The lectures on the control theory", 1975. p. 60. 
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2 Domains of asymptotic stability 



Let us consider the differential system 

x= f(x), (1) 
where x = (x±, X2, x n ) G W 1 and 

f(x) = (f 1 (x)j 2 ( x ),...j n (x)y. 

The vector-polynomials f p (.) and their coefficients a P h 

,h,...,l n sa tisfy the conditions in the 

Introduction. 

We assume that f(x) ^ for all x ^ in some neighborhood of the origin of coordinates 

0. 

We will call the order deg{f p {.)) of the polynomial f p (-) the maximal degree of the 
polynomial f p {.) in the variables Xj,j £ 1 : n, in totality i.e. 

deg(f p (x)) = max (h + l 2 + ... + l n ). 
So if for n = 2 and kj / 

fl(x) = l\\x\ + lvix\ + ^13^1^2, f2(x) = hixf + 122X2 + I23X3, 

then the order of fi(x) is equal to two and the order of f2(x) is equal to three. 

We will call the order of the function /(.) the maximal degree of the polynomials 
fp( x ),p £l:n, regarding the variables Xj,j £ 1 : n, i.e. 

deg(f(x)) = max deg(f p (x)). 

Consequently, the order of the function /(.) for the example written above is equal to 
max(2, 3) = 3. 

We assume that the definitions of stability and asymptotic stability are known (see 
[l]-[4])- 

Definition 2.1 [1], [2] A domain D, G intD, consisting of whole trajectories of the 
system (1), is called a region of asymptotic stability if the limit 

|| x(t,x ,t ) (2) 

when t — > 00 is fulfilled for any initial point xq € D and any solution x(-),x(to) = xo, of 
the system (1). In this case we will call the system (1) asymptotic stable in D. 

The following problem is important for practical applications. 
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Problem 2.1 It is known that zero solution of the system 



(3) 



A[n xn],x £ W 1 , is asymptotic stable if all eigen-values of the matrix A have negative real 
parts. In this case the system (3) is asymptotic stable in 3? n . 
Consider now the differential system 



where ip(-) is a vector polynomial with degree not less than two, <^(0) = 0. 

What are the conditions on the coefficients of the vector- polynomial </?(•) from (4) 
under which the system (4) is asymptotic stable in a given domain 0,0 £ intD, i.e. for 
any solution of the system (4) and any initial point xo G D the limit (2) is true? 

The considered circle of questions includes the problem about center and focus. This 
problem is formulated in the following way. 

Problem 2.2 Assume that the origin of coordinates (0,0) for the system (3), n = 2, 
x G %t 2 is the center. The necessary and sufficient conditions for this are that all eigen- 
values of the matrix A are imaginary. By adding a vector polynomial </?(•) (system (4)) 
the center (0,0) can be a focus. It is needed to find conditions on the coefficients of the 
polynomial </?(•) that the point (0,0) was the center of the differential system (4). 

Let us rewrite the system (1) in equivalent form 



The elements aij(x) of a matrix A(x)[n x n] are continuous polynomial functions of x. 
Conversion (4) to (5) is not unique. It can be done in an infinite number of ways. In fact, 



x= Ax + tp(x) 



(4) 



x= A(x)x. 



(5) 



Oijix) = 



E 



a 




n 



{x)x l ix l 2...x) 



n 



■n 



h+h+-+l n <deg(f) 



if lj > 1. We have the following correlation for the coefficients 




(6) 
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for all x from some region D, € intD. 

For instance, let the system (1) have the right part 




Then 




a\xix 2 + ct2x\ 1 + j3ix\ + /3 2 2qx 2 

+ 72^2 8\X\ + 82X1X2 
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where the coefficients ccj, (3i,ji, = 1,2 such that 

ol\ = 1, 71 + <5i = 3, 

a 2 + (3 2 = 2, 72 + <5 2 = -2. 
The system (4) for the given vector function /(•) can be rewritten in the following form 

x= A x + C(x)x (7) 

where 

,2 R.J2 



r<( \ — ( ol\X\X2 + a 2 Xr; Pixf + P2X1X2 

1 71^1^2 + 72^2 <5l x l + ^2XiX 2 




A = 

Since the eigen-values of the matrix Aq are X\ t 2 = ±i, the point (0,0) is the center for the 
linearized system. But it is question for the nonlinear system (7). 

We will consider all possible continuous matrices A{x) whose elements are polynomial 
functions of x for which the system (5) is equivalent to the system (1) . We will denote 
the set of all such matrices by A. 

Let us solve the problem 2.1. 

Theorem 2.1 In order that a domain D consisting of the whole trajectories of the system 
(1) i.e. x(-,Xo,to) £ D,xq € D, for all t > t^ is a region of asymptotic stability it is 
necessary and sufficient that there is a matrix A(-) £ A of the system (5) in the domain 
D whose eigen-values have negative real parts at any point x G D, x 7^ 0. 

Proof. Necessity. Let a domain D be asymptotic stable. Consider any trajectory 
xq, to), x(to) = xq. There is such a transformation £ = X{x) of W l in a neighborhood 
of a point x\ = x(t\) that transforms the system (1) to the differential system 

k= B{m (8) 

where the matrix B has the eigen-values with negative real parts. 

Indeed for any point x(t\, Xo, to),t\ > to there is a linear transformation T x defined in a 
small neighborhood of the point x(ti, xo, to) and that is the linear part of the transforma- 
tion X(.) in this neighborhood, so that the system (1) is transformed by this transformation 
to the differential system (8) and the vector 

xi =x (ti,x ,to) 

is transformed to the vector 

£1 =£ (ti,x ,x ). 

Any differential equation can be defined by its current of tangent vectors. If we choose an 
asymptotic stable linear differential system like the system (3) which current of tangent 
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vectors is close to current of tangent vectors of our system locally for each point x of the 
domain D then the eigen-values of the matrix of this asymptotic stable linear differential 
system will have negative real parts. Consequently, the above mentioned transformation 
X{x) exists. 

As soon as the system (8) is transformed to the system 

x= T^B^TxX 

under the linear transformation 

i = T x x 

in a neighborhood of the point £1 and all locally linear transformations of the system (1) 
have the form (5), then the matrix A(-) G A exists for which 

A(x(t 1 ,x ))=T- 1 B(Z 1 )T x (9) 

and the eigen-values of the matrix A(x(t±)) are equal to the eigen-values of the matrix B. 

Since the transformations T x and T~ x are continuous with respect to x, the matrix 
A(.) is also continuous with respect to x G D . 

Sufficiency. Let the system (1) admit such a transformation to the form (5) that all 
eigen-values of a matrix A(x) have negative real parts for all x G D, x ^ 0. We prove that 
the domain D is a region of asymptotic stability. 

Fist we remark that in general case not any solution of the system (1) 

x(t) = x(t, x , to), x(t ) = x 

can be represented in the form 

x (t) = e Jt « x (10) 

where integration is taken along the integral curve x(.). 

Instead of (10) we will use the sequential approximations {xk{t)} of the system (1) on 
the segment [to , t] having the form 

x k (t) = e A ^-^ t ~ t ^x k . 1 (t) (11) 

where {U} is a subdivision set of the segment [to,t] and 

x i+l = e A{x ^ u+1 - u) Xi, i€0:(fc-l). (12) 

If we substitute (12) into (11) we will get 

Xk (t) = e A ( x k-i)(t-tk-i) e M x k-2)(tk-i~tk-2) ...e^oX*! - *")^. (13) 

Since A(xi) — > x _^ x ^ A(x(t)), where x(.) is a solution of (5), we have 

U -> T 



(14) 
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uniformly on r £ [to, t] when k — > oo and 

maxj G i : fc | ti — ti-i \— >k 0. 

It is obvious that all matrices A(xi),i G : (k — 1), have the eigen- values with negative 
real-valued parts for enough big k. 

Denote by \i(xj) the eigen-values of the matrix A(xj),i G 1 : n,j G : k — 1. Then a 
number C > exists that 

|| e A(a; fc _ 1 )(t-t A ._ 1 ) e A(a; fc _2)(t fc _i-t fc _2) _ _ e A(xo)(ti-to) Xo || < (J e Ht)(t-t ) || 

where 

A(t) = lim max i?e Aj(xj). 

fe^oo j e o : fe - 1, 

i G 1 : n 

Since A(i) < for any t > to, 

X(t)(t-t ) < 0. 

Let 

-a < X(t)(t-t ) < 0,a > 0. 

This relation means that the trajectory x(-) can not go to a stationary point or be a 
stationary orbit. Indeed if the trajectory x(.) has a stationary point or is a stationary 
orbit then 

lim X(t)(t — to) = — oo. 

t—>co 

Consequently, 

|| x(t,x ,t ) || — > 

when t — > oo. We obtain the contradiction. Thus, the origin of coordinates can be only 
a stationary point. The sufficiency and the theorem are proved. □ 

Remark 2.1 It follows from the Necessity of the Theorem 2.1 that there is one-to- one 
correspondence between the matrices A(x) G A of the system (5) and locally linear trans- 
formations T x of the system (1) defined in a neighborhood of the point x(see (9)) . Indeed, 
there is a matrix A{.) for any transformation T x and, on the contrary, there is own coor- 
dinate system for any matrix A{x) G A, in which the equation (9) is true. 

Corollary 2.1 .It follows from the Necessity of the Theorem 2.1 that the degree on x of 
any element of a matrix A(x) which we denote by aij{x),j G 1 : n, that is a polynomial of 
x, does not exceed the degree of the vector-function /(•). 

Proof The sum of the coefficients a,ij(x) ■ Xj over j is equal to zero for those elements 
that do not belong to /«(•)• And on the contrary the sum of the coefficients ^ ln 
over j is not zero for those elements that are the terms of fi{ ). It follows from here that 
choosing ^ ^ from the system (6) we can only consider the elements of the vector- 
polynomial /(•). □ 
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Corollary 2.2 . For asymptotic stability of a solution x(.,xo,to), x(to) = xq, of the 
differential equation (1) it is sufficient that there was 5(to,e) > for any to,s > that for 
II x o ||< 8 



Proof. As soon as the eigen-values of a matrix A(-) have negative real parts in any 
neighborhood of the origin of coordinates from which any solution tends to zero-vector , 
then for any k we have from (13) that 



i.e. all solutions are in an e -neighborhood of the origin of coordinates. The latter means 
stability. The Corollary is proved. □ 



3 To problem about center and focus 



The following problem is interesting from technical point of view: to recognize focus or 
center i.e. to give the conditions when trajectories turn around and go to the origin of 
coordinates or remain closed. 

From the start we will give the definitions of focus or center. Beforehand we define 
trajectories turning around some ray with the initial point € !R n infinitely often. 

Definition 3.1 Let us say that a trajectory x(t,xo,to) turns around a ray I G 3?" with the 
initial point € ffl 1 infinitely often if the radius vector r(t) = x(t,XQ,to) forms with the 
ray I an angle ip(t) taking all values from a segment [991,992] infinite number of times. 

Definition 3.2 We will call the point e focus of the system (1) if 

1 ) the point is asymptotic stable ; 

2) any solution x(i,xo,io) of the system (1) turns around a ray I £ ffl 1 with the initial 
point OgSR" infinite number of times. 

Definition 3.3 We will call the point a center of the system (1) if 

1 ) the point is stable ; 

2) all solutions x(t, xo,to) of the system (1) remain in some neighborhood of the point 
and are closed loops. 





Xk(t) \\<\\ x ||< e 



and in limit on k 



x(t,x ,t ) \\<\\ xq \\< e 
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Definition 3.4 We will call the point G K ra a node of the system (1) if 

1 ) the point is asymptotic stable ; 

2) there are not a ray I G ffi 1 with the initial point G W l and a solution x(t,Xo,to) of 
the system (1) turning around the ray I infinite number of times. 

Theorem 3.1 In order that the point £ S" is a focus of the system (1) with a right 
polynomial part /(.), f(x) / for x ^ it is necessary and sufficient that 

1. There is such a continuous matrix A(-) £ A of the system (5) that all its eigen-values 
at any point x,x ^ 0, from some neighborhood D of the origin of coordinates, con- 
sisting from the whole trajectories, have negative real parts and non-zero imaginary 
parts. 

2. There is no matrix A(-) G A of the system (5) with the negative real-valued eigen- 
values at all points x G D, x ^ 0. 

Proof. Necessity. Let the point be a focus for the system (1). There is such a 
transformation £ = X(x) so that the system (1) can be rewritten in a neighborhood of 
some point x\,x\ / 0, in the form 

where B(£) is a matrix whose eigen-values have negative real parts and nonzero imaginary 
parts. 

Consider any trajectory x(-, xq, to), x(to) = x$. Then for any point X\ G x(t, x$, to), x\ ^ 
0, there is a transformation T Xl in a neighborhood of the point X\ that the system (1) can 
be rewritten in the form indicated above. 

If the trajectories of the system (1) do not turn around any ray I G K n with the initial 
point G 5ft™ infinite number of times then such system can be transformed by some 
continuous transformation that corresponds to some matrix A(-) G A that is the matrix 
of the system (5) whose eigen-values are negative real- valued numbers at all points x G D 
(theorem 2.1). The latter is impossible according to the Condition 2 of the theorem. 

Sufficiency. Let the conditions 1,2 of the theorem be true. Prove that the point is 
a focus of the system (1). 

Any solution x(-, Xo, to), x(to) = xq, is the limit of some sequence x^{t) on k — > oo 
obtained from (13). According to the conditions all matrices A(xj),j G 1 : k, have the 
eigen-values \\{xj),j G 1 : k,l G 1 : n with negative real- valued parts a/(-) and nonzero 
imaginary parts bi(-) : Xi(x) = a\(x) + ibi(x),i 2 = —1, I G 1 : n, that are nonzero at any 
point x G D. It follows from here that the point is either focus or node. 

We will prove that under the Condition 2 all trajectories turn around a ray I G K n 
with the initial point G 3?" infinite number of times. If it is not true, then there is such 
a transformation which corresponds to a matrix A(-) G A whose eigen-values are negative 
real-valued numbers for x G D, x / 0. This contradicts the Condition 2. The sufficiency 
and the theorem are proved. □ 
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Corollary 3.1 If instead of the Condition 2 of the Theorem 3.1 we require that there is a 
matrix A(.) £ A of the system (5) with the eigen-values \%{x), i G 1 : n, described in the 
condition 1 for which 

Im Xi(x) > a, a > 

or 

Im \i(x) < a, a < 0, 

where imA&(.) denotes the imaginary part of Xk(-), for x G D,x ^ 0, then all trajectories 
turn around a ray I G 3ft ra with the initial point G 3ft ra infinite number of times. 

Proof. From the representation of any solution of the system (1) as a limit of sequence 
{xj(.)} from (13) it follows that a solution x(-,xo,to) has a finite number of turns around 
any ray I G W 1 with the initial point G W 1 if and only if for any i G 1 : n the sum 

53 ImXi(xj-i)(tj -tj-i), 
jei-.k 

where Xj(.) was defined in (13), or in the limit on k — > oo the integral 

/ Im Aj(x(r))dr 

is convergent. But under the condition of Corollary 3.1 this integral is divergent. The 
corollary is proved. □ 

Mathematicians were concerned about the problem of recognition of center ore focus for 
a long time. The theorem given below for two dimensional system with a right polynomial 
part has the point = (0, 0) as a center. 

Theorem 3.2 In order that the point = (0, 0) be a center for the two-dimensional 
system (1) it is necessary and sufficient that there was a matrix A(.) G A of the system 
(5) whose eigen-values are non-zero imaginary numbers in a neighborhood S of the origin 
of coordinates, consisting from whole trajectories, where f(x) ^ for x ^ 0,x G S. 

Proof. Necessity. Let the point = (0, 0) be a center of the system (1). Prove that the 
condition of the theorem is true. 

Take such a small neighborhood S of the origin of coordinates where all trajectories 
starting in S are closed loops. Consider any trajectory x(-, Xo, to), x(to) = xq. Take a point 

x 1 = x(ti,x ,t ). 

There is such a transformation T x in a small enough neighborhood of the point x±,x± ^ 0, 
that the system (1) can be rewritten in the form 

and the matrix B(.) has the imaginary eigen-values for all £ G S. As was proved above 
(Theorem 2.1, the Proof of the Necessity) the matrix B(-) corresponds to some matrix 
A(-) that like the matrix £>(•) has the imaginary eigen-values. The necessity is proved. 
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Sufficiency. Let the condition of the theorem be true. Prove that the point = (0, 0) 
is the center. 

As soon as any transformation T x is defined in a small neighborhood of the point x, 
then the above statement is not sufficient for the point = (0, 0) to be a center. The point 
= (0, 0) may happen to be a focus convergent (nonconvergent) with a finite number of 
turns around any ray I G 5ft n with the initial point 0. But it is impossible because from 
the limit of the sequence (13) we conclude that for any k 

II x k(t) || = || x || • 

From (14) it follows that the solution x(-,xo,to) can't go to the origin of coordinates. 

A set of closed and not closed loops could occur. From the representation of an integral 
curve x(-,xo,to) in the form (13) and (14) it follows that there is not a matrix A(-) G A 
with the imaginary eigen-values at all points of a small neighborhood of the origin of 
coordinates. The latter contradicts the condition of the theorem. 

There is another case when for all i G 1 : n the sum 

ImXiixj-^itj - tj-_i) 

jei-.k 

or in the limit on k — > oo the integral 

/ Im Xi(x(r))dT 
J t 

has a finite value where Aj(-), i G 1 : n, are the eigen-values of the matrix A(-). In this 
case the trajectory x(-,xo,to) goes to some stationary point x G S, x ^ 0. It means that 
f(x) = 0. This contradicts the condition of the theorem. Consequently, the point is the 
center. The sufficiency and the theorem are proved. □ 

Corollary 3.2 If a matrix A(.) G A of the system (5) exists with the imaginary eigen- 
values \i(x) for all x G S where S is a neighborhood of the origin of coordinates 0, 
consisting from whole trajectories, and 

Im Xi(x) > 

or 

Im Xi(x) < 

for all x G 5, x ^ 0, and i G 1 : n then the point is the center of the n- dimensional 
system (1). 

Proof. Indeed, if the conditions of the Corollary 3.2 hold, then the conditions of the 
Theorem 3.2 hold as well. □ 

The Corollary 3.2 does not require that f{x) ^ for all x ^ from some neighborhood 
of the origin of coordinates. 
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Remark 3.1 The theorem can be proved for any n- even- dimensional spaces. It is not 
difficult to see that for n dimensional spaces with n odd, n = 2k + 1, k is a natural number, 
the system (1) can not have the origin of coordinates as a center. Indeed, there is no 
imaginary number among the eigen-values of any matrix A(-) for odd n. 

Remark 3.2 Let us consider the differential system 

x=f(x,t) (15) 

where x = (x±, X2, x n ) G K n and 

f(x,t) = (f 1 (x,t),f 2 (x,t),...,f n (x,t))*. 

The vector- function f p (.) is a polynomial of x i.e. 

fp( x ) = X! a J^2,.--,*n(*) X l lx 2 — x n > 

p G 1 : n, li, I2, ■■■l n are non-negative integers, a^\ 2 1 (t) are continuous real-valued func- 
tions and * is the transposition sign, I p is a finite set of indexes of the polynomial f p {.). We 
will assume that f(x, t) ^ for all x ^ in some neighborhood of the origin of coordinates 
and t > to. 

In this case we will denote by \i(x,t),i G 1 : n, the eigen-values of a matrix A(.). If 
we demand that all cited above statements about the eigen-values Xi(x,t),i G 1 : n, are 
true for all t > to then we obtain the sufficient conditions for the cited above theorems and 
corollaries. 

Problem 3.1 For the system 

x=f(x,t) (16) 

where /(., .) is a polynomial of x and t it is required to find some conditions when the 
system (16) is asymptotic stable in a domain D,0 G intD. 

The ideas stated above do not apply to the system (16) because the terms of /(., .) 
can be unbounded along some solution x(-,xo,to) in D. If we require the terms of /(., .) 
to be bounded along the trajectories of the system (16) and that \\x (i,xo,io) II goes to 
zero uniformly on xq when t — > 00, then this system can be transformed to an equivalent 
stationary system. This idea will be developed in following articles. 



4 Stabilizing control 



Let us go to the question of finding a stabilizing control. 
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Let us consider the differential system 

x=f(x,u) (17) 

where x = (x±, X2, x n ) G ffi 1 is a phase vector, u = (ui, u 2 , u r ) G ffi is a control, 
f(x,u) = (/i(x, u), f 2 (x, u), f n (x, u))* is a vector-polinimial of x and u with constant 
real-valued coefficients, i.e. 

fp(x,u) = ^ a h!l2,...,l n ,m 1 ,m 2 ,...,m r X l X 2 — x n U™ 1 U™ 2 ■ ■ ■U™ 7 ' , 

h,h,---,ln,mi,m2,...,m r £l p 

(p) 

p G 1 : n, l 1 ,l 2 ,...,l n ,m 1 ,m 2 ,...,m r are non-negative integers and af^ ^^^ mr 
are real- valued numbers, I p is a finite set of indexes of the polynomial f p (.). 

Let us assume that the zero-vector = (0,0,..., 0) G ffl 1 is a solution of the system 
(17) for u = G 5? r . 

Definition 4.1 The control u{x) = (u±(x), u 2 {x), u r {x)) G 5ft r is called stabilizing 
in a domain D G K n , G intD, for the system (17) if any solution of (17) x{t) = 
x(t, xo, to, u(x(t))) satisfies the limit 

|| x(t,x Q ,t ,u) || — > (18) 

when t — > oo, x(to) = xo £ D. 

As stated above the condition (18) is sufficient for the system (17) to be asymptotic stable 
in D, i.e. the zero solution of the system (17) is asymptotic stable and the limit (18) is 
true for any initial point x(io) = xq £ D. 

Problem 4.1 It is required to find such a stabilizing control u = u{x) in a given domain 
D,0 G intD,x G D that for any solution x(t,xo,to,u),x(to) = xo of the system (17) the 
correlation (18) was true. 

For linear system 

x= Ax + Bu, (19) 

where B[n x r] is a matrix of amplification coefficients, a stabilizing control u G W can 
be chosen in the form u = Cx so that (18) is true. 

Theorem 4.1 [1],[4]- If the rank of the system 

B, AB, A 2 B, A n ~ 1 B 

is equal to n then we can always construct a stabilizing control in !R n in the form 

u = Cx, 

where C[r x n] is some matrix. 
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Since the matrix B is denned by technical essence of the system, we can choose it con- 
structing a system by ourselves . Therefore the condition of the theorem can be always 
satisfied. 

The theorem was written for analogy and comparison with the results below. 

Our goal is to construct a stabilizing control u(x) that a domain D were a domain of 
asymptotic stability. In common case it is not always can be done. But if we change the 
system (17) a little bit the problem can be solved. Instead of the equation (17) we will 
consider the equation 

x= f(x,u) + <p{u) (20) 

where </?(•) is a vector-polynomial (f(u) = ((pi(u),ip2(u), ...,(p n (u))* with the degree not 
bigger than the degree of the function f{z) as a function of z = (x,u), 

ii,i2, ■■■iveAfp 

(p) 

where p G 1 : n and of l i2 ir are constant real- valued numbers, i±, 12, i r are non- negative 
integers, M p is a finite set of indexes of the polynomial <f p (-)- 

In practice it is possible to construct </?(•) because we choose a control u(-) ourselves. 

We will look for a stabilizing control u = u(x) in form of polynomial in x. 

Theorem 4.2 For any domain D,0 £ intD, the vector-polynomials u(.) and </?(.) can be 
chosen such that 

1) D is a region of asymptotic stability for the differential system (20); 

2) the degree of u(x) does not exceed the degree of the vector-polynomial f(x,u) as a 
function of x; 

3) the degree of <£>(■) ^ s n °t bigger than the degree /(.) as a function of 
z = (x, u). 

Proof. We will prove that a vector polynomial u{x) can always be chosen satisfying the 
following conditions 

1) the degree of u(x) does not exceed the degree of the vector-polynomial f(x,u) as a 
function of x; 

2) u(x) is a stabilizing control for (20) in the domain D. 

We will use the results obtained before (Theorem 2.1). 

Let us substitute in (20) the control u = u(x) in the form of a vector- function of x. 
The system (20) is rewritten as 

x= f{x) (21) 

or 

x= A(x)x. 

We can write conditions for the matrix A{x) to have the eigen-values with negative real 
parts. 

There is no difficulty calculating the degree of the function f(x,u) and the function 
<p(u) as a function of x after substituting u = u(x). 
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Let us denote by l x and l u the degrees of the function f(x,u) in the variables x and u 
correspondently. Then after substituting u = u{x) the degree of the function /(., .) as a 
function of x is not bigger than l x + l u l x and the degree of the function (p(.) as a function 
of x is not bigger than l x (l x + l u ). It is easy to see that 

lx ~\~ Ixlu ^ Ixifx ~\~ lu)- 

Having chosen the function <p(-), the coefficients of the vector polynomial f(x) can be 
chosen so that the matrix A(x) has the eigen-values with negative real parts. The theorem 
is proved. □ 

In an analogous way it can be proved the following theorem. 

Theorem 4.3 For the differential system (20) with even n the vector-polynomials u(-) 
and f(-) can be chosen such that 

1) the origin of coordinates is the center for the system (20); 

2) the degree of u(x) does not exceed the degree f(x,u) as a function of x 

3) the degree of tp(.) is not bigger than the degree of /(.) as a function of z = (x,u). 

This theorem can be used in physics of plasma specially for stabilization of plasma in 
reactor. 

(v) 

Remark 4.1 In the case when the coefficients t m2 mr depend on t the coef- 

(P) !))>!!! 

ficients of l i2 ir will depend on t as well. We can try to find these coefficients that the 
conditions of the theorems of the previous section were true for all t > to- 



5 One aspect of application 



We will study for instance how to choose a control u(-) so that a given domain D consisting 
from the whole trajectories of a differential system were a region of asymptotic stability. 
For that we have to solve the following optimization problem. 

Let us substitute a vector-control u(-) in a form of polynomial of x with a degree m 
not bigger than the degree of the the vector-polynomial f(x,u) as a function of x in the 
equation (20) and rewrite (20) in the form (21). Denote the eigen-values of a matrix 
A(-) in (21) by Aj(s,a ij -,dj[ ) )i2) ..., in ,&S 2 , ...,i r ) where a H and d h,i 2 ,...i n are the coefficients 
of the matrix A(-) and the vector- function u(-) correspondingly, b!?\ i , p G 1 : n, are 
the coefficients of the vector- function </>(.). Then our problem is reduced to the following 
optimization problem: to find out such continuous functions otij(x),i,j € 1 : n, x £ D and 
the numbers , b^\„ ,- , p € 1 : n, that the following correlation 

A,i,-.o, ;! ,-).<;, : , ;( ). (22) 
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was true for all x G D,x ^ 0, where aij(x),i,j G 1 : n, are connected with each other by 
linear equations for each x <G D, x / 0, by other words we should solve the next problem 



m inf SU P i c nf ReX j( x ^ a ij^A,i2,...,i^K,i2,...,ir) < V(5 > 0, 

ft riir ,P6l:n x€D\B?(0) 



where Bg(0) = {z£ 1 1[ z \\< 5}, S— is any sufficient small number for which D D Bg(0). 

It is obvious that the inequalities (22) can be replaced by an equivalent system of 
inequalities for coefficients of the characteristic polynomial of the matrix A(.). To solve 
this system it is easier than to find the eigen- values of the marix A(.). 
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